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QCD undergoes a deconfining transition at high temperature to a "quark-gluon plasma" phase 
where hadrons may become unbound. In this work, meson correlation functions at non-zero mo- 
mentum are studied both in the confined and plasma phases using the Maximum Entropy Method. 
In particular, both the longitudinal and transverse modes of the vector correlation functions are 
considered. Only in the case of light quarks in the plasma phase, we find that both longitudinal 
and transverse spectral functions have a non-zero intercept at zero energy. 
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1. Introduction 

There has been considerable interest in the quark-gluon plasma phase of QCD over recent 
years. On the experimental side, experiments at RHIC have led to an increased understanding of 
this phase of matter. This has been matched by a similarly intense theoretical understanding using 
e.g. hydrodynamic models and lattice simulations. Despite this considerable work, there remains a 
great deal of uncertainty over many properties of the quark-gluon plasma. In particular, the Particle 
Data Book does not contain a single entry for this phase [|l|]. 

This work aims to extend our knowledge of the physics of QCD in the plasma phase by study- 
ing the properties of mesonic spectral functions, p(co), both below and above the deconfining 
temperature, T c . The zero energy limit, ft) — > 0, of spectral functions gives information on hydrody- 
namic structure and transport coefficients; we will be particularly interested in the conductivity, a, 
and the related diffusivity, D = c/S, where S is the charge susceptibility. These can be obtained 
from vector spectral functions. In this paper, we outline an extension of our earlier work [§] to the 
case of non-zero momentum |]3|l. This allows us to study the longitudinal and transverse compo- 
nents in the vector channel which can, in principle, unlock interesting information about transport 
and hydrodynamics [0]. For more details on transport coefficients and lattice QCD, we refer to the 
reviews [§,§]. 

2. Lattice and Fitting Details 



We begin by defining the usual spectral function in terms of the Euclidean two-point function, 

G(t,p), 

G(t,p) = J^p(co,p)K(t,co), 
where p denotes the momentum. The kernel is given by 

cosh[fl)(f - l/(2T))} 



K{t,co) 



sinh[fti/(2r)] 



where T = \/{aN t ). We note that the extraction of a spectral density from a lattice correlator is 
an ill-posed problem, since the correlator, G(t), is known at only ^(10) time points, whereas the 
spectral function, p(ft)), is, in principle, a continuous function. The usual solution of this problem 
implements the Bayesian analysis technique of the Maximum Entropy Method (MEM) [{7|]. 

At finite temperature the kernel is singular and independent of euclidean time in the limit that 
ftw0[§: 

lim K((D,t) = — + 0((o). 
eo-S'O ft) 

In Ref. [Q] we uncovered that this singularity affects the reliability of the MEM procedure near 

(O = 0. Fortunately it can be trivially corrected by defining a rescaled kernel and spectral function, 

K((o,t) = —K{(o,t), p(ft>) = — p(ft>), 

and performing MEM on 

G(t) = f d ^p((o)K((0,t). 
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Cold Hot 



Spatial Volume 


N^ x N t 




48 3 x 24 




64 3 x 24 


Lattice spacings 


a- 1 




~4GeV 




~ 10 GeV 


T 


l/(aN t ) 


T ~ 160MeV~ 0.62T C 


T r 


-420MeV~ 1.5r c 


Statistics 


N cfg 




-100 




~100 


Quark Masses 


ma 




0.01 & 0.05 




0.01 & 0.05 



Table 1: Lattice parameters used in the simulation. Estimates for the lattice spacing and temperature are 
taken from Ref. [|I§. 

We demonstrated that this redefinition removes the unphysical behaviour of the spectral functions 
constiucted with MEM near CO = 0, present when the tiaditional kernel is employed [g]. 
In the MEM procedure, we used the default model 

m(co) = rao(l + aco) 

with mo determined from a fit of Jdcom(co)K(co) to the data. This definition of Tn(co) matches 
the expected perturbative behaviour p((o) ~ cop(co) ~ CO 2 as CO — > °° and also allows a non-zero 
intercept in p/co ~ p (corresponding to a transport peak) as CO — > 0. For a study of lattice spectral 
functions at nonzero momentum in the infinite temperature limit, see Ref. [^]. 

The lattice action and parameters used are identical to those in Ref. [g], i.e. a simple Wilson 
plaquette action with quenched, staggered fermions. The lattice parameters used are displayed in 
Table |l|. In Ref. [Q] we determined the electrical conductivity, a/T = 0.4 ±0.1, a result which 



has since been confirmed [11]. We extend the work in Ref. [g] further by using twisted boundary 
conditions with the same choices of twists as in Ref. []l2|]. This enables us to access a large range 
of momenta which are listed in Table ||[ Preliminary results can be found in Ref. (J3J]. 

Since we use staggered correlators which have contributions from both parity partners, 

dco 



G(t) = J ^K(t,co)[p(<o)-(-l) t/a p((o) 



we are forced to apply MEM to even and odd times separately, and then reconstruct the physical 
spectral function from 

3. Longitudinal and Transverse Correlators 

In Fig. |] we show the vector correlators for both the longitudinal and transverse modes, Gl,t. 
for momenta pL = 2, n, 27T, 3n as a function of euclidean time. We display these as a ratio with the 
average correlation, GAve = (Gl + 2Gt)/3. Note that we use local (rather than conserved) currents 
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pL 


\P\L 




Longitudinal 


Transverse 


(0, 0, 0) 

(2, 0, 0) 

(0, ±7T, 0) 




2 

K 


3.72 
>8 

^6.59 
12 
^9.63 


v 2 
v 2 
v 2 


V 2 & V 3 

Vi & v 3 


(-2, ±7T, 0) 
(0, ±27T, 0) 


27r~6.: 


v 3 

Vi & v 3 


(2, ±27T, 0) 
(0, 37T, 0) 


2Vl±^ 2 - 

37T~9.' 


v 3 
V\ & v 3 


(-2, 3tt, 0) 


V4±9tt 2 - 


v 3 



Table 2: A list of the momenta studied showing which components of the vector current, V, = y/YiV' are 
longitudinal or transverse. The momenta components pL which are multiples of n are obtained from the 
usual Fourier sum, and those pL components which are integer-valued are obtained from twisted boundary 
conditions. 



a 






a 



1.4 
1.3 
1.2 - 
1.1 
1 
0.9 
0.8 - 



i 1 1 r 



i 1 1 r 
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Figure 1: Longitudinal and transverse vector correlation functions normalised by the average correlation 
function, Gavs = (Gl + 2Gt)/3, for various momenta, p, as a function of euclidean time. Data points for 
each momentum (for a given time) are offset horizontally for clarity; from left to right they are \p\L = 2, n, 
2% and 3k. 
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Figure 2: Longitudinal (left) and transverse (right) vector spectral function predictions from Ref. 
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Figure 3: Longitudinal (left) and transverse (right) vector spectral function for the light (ma = 0.01), hot 
case. The main plot shows the CO ~ region, and the insert the full CO range. 



in our correlators definitions. As can be seen, there is a clear difference between these modes; the 
longitudinal correlator is consistently larger than the transverse for each time slice and momentum 
value. It would be interesting to understand this feature analytically. Note that the momentum 
dependence is much stronger for even timeslices. This is an artefact of the staggered formulation. 

The diffusivity, D, can in principle be obtained from the momentum dependence of the vector 
spectral function (at small mass). In Ref. [Q] a prediction was made of the energy dependence of the 
vector spectral function for both the longitudinal and transverse modes (for the massless case in the 
plasma phase). This is shown in Fig. ||. As can be seen, their prediction for the co — > behaviour 
differs between these two modes: Pl(co = 0, p) = for p / whereas pi{co = 0,p) / 0. 

In Fig. ||, we show the spectral functions obtained via MEM for the longitudinal and transverse 
vector case for quark mass ma = 0.01. As can be seen, there is a clear non-zero intercept for both the 
longitudinal (contradicting [Q]) and transverse modes. The spectral function for the heavier quark 
mass, ma = 0.05, is shown in Fig. |J We note that the intercept is zero for both the longitudinal 
and transverse cases for this heavier quark mass yielding a null result for the heavy-quark diffusion 



coefficient [14]. 
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Figure 4: Longitudinal (left) and transverse (right) vector spectral function for the heavy, hot case. The 
main plot shows the co ~ region, and the insert the full CO range. 




Figure 5: Pseudoscalar spectral function for the ma = 0.01 (left) and ma = 0.05 (right) cases. The main plot 
shows the CO ~ region, and the insert the full CO range. 

We checked the spectral function for both the quark masses ma = 0.01 and 0.05 in the cold 
phase, and found a zero intercept for both longitudinal and transverse modes. This agrees with our 
expectation that there are no transport features in this phase. 

In Fig. Hi the pseudoscalar spectral functions are shown in the hot phase. We observe a zero 
intercept for this channel again in agreement with the expectation that there are no transport fea- 
tures. 

4. Conclusions 



In this work we have extended our previous study |Q] of the spectral functions of vector corre- 
lators in the cold and hot phase of QCD to non-zero momenta using twisted boundary conditions 
to allow a finer momentum resolution. We are thus able to calculate and contrast the longitudinal 
and transverse modes of the vector two-point function. We find that a zero intercept at zero energy 
for the spectral function in all cases except for the hot, light quark vector mesons where we find 
p (co — > 0, p) 7^ for both the longitudinal and transverse modes. This is despite the fact that the two 
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modes' two-point correlation functions have a different qualitative behaviour. Our future plans are 
to study the momentum dependence of these modes with the aim of further studying hydrodynamic 
structure at zero and nonzero momentum [|l3|]. 
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